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We report on a systematic study of two dimensional, periodic, frustrated Ising models with a quan- 
tum dynamics introduced via a transverse magnetic field. The systems studied are the triangular and 
kagome lattice antiferromagnets, fully frustrated models on the square and hexagonal (honeycomb) 
lattices, a planar analog of the pyrochlore antiferromagnet, a pentagonal lattice antiferromagnet as 
well as a two quasi one-dimensional lattices that have considerable pedagogical value. All of these 
exhibit a macroscopic degeneracy at T = in the absence of the transverse field, which enters as 
a singular perturbation. We analyze these systems with a combination of a variational method at 
weak fields, a perturbative Landau- Ginzburg- Wilson (LGW) approach from large fields as well as 
quantum Monte Carlo simulations utilizing a cluster algorithm. Our results include instances of 
quantum order arising from classical criticality (triangular lattice) or classical disorder (pentagonal 
and probably hexagonal) as well as notable instances of quantum disorder arising from classical 
disorder (kagome). We also discuss the effect of a finite temperature, as well as the interplay be- 
tween longitudinal and transverse fields — in the kagome problem the latter gives rise to a non-trivial 
phase diagram with bond-ordered and bond-critical phases in addition to the disordered phase. We 
also note connections to quantum dimer models and thereby to the physics of Heisenberg antifer- 
romagnets in short-ranged resonating valence bond phases that have been invoked in discussions of 
high-temperature superconductivity. 
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I. INTRODUCTION 

The study of frustrated magnetic systems began a 
half centugy ago with the realization by WannierEI and 
HoutappeB that the antiferromagnetic Ising model on a 
two-dimensional triangular lattice does not order down 
to the lowest temperatures and exhibits a finite entropy 
density even at T — 0, in contrast to the naive expecta- 
tion from the third law of thermodynamics. These two 
observations are related, and can be traced to the frus- 
trated nature of the couplings even at the level of a single 
plaquette (Fig. |l|) where one sees that it is not possible 
to minimize the energy of all three bonds simultaneously 
leading to six (instead of the two for ferromagnets) min- 
imum energy configurations. From this we may infer a 
macroscopic entropy density for the triangular lattice at 
T = and rationalize the absence of any ordering by the 
easy accessibility n of a large number of configuration at 
all temperatures .S 

The twin observations, of a non- vanishing entropy and 
lack of order, when contrasted with the ordering tran- 
sition and small number (two) of ground states in the 
ferromagnetic Ising model on the same lattice, typify the 
striking behavior of such "maximally" frustrated classi- 
cal models whose catalog, by now, includes also models 
with continous spins such as the Heisen.berg magnet on 
the three dimensional pyrochlore lattice.tffl Other, by def- 
inition, less, frustrated models exhibit non-obvious phase 
transitions^ at finite temperatures. These may involve 
a singular privileging of the ground state manifold as a 
whole |-as in the case of the XY magnet on the kagome 
lattices or, more typically, the elegant phenomenon of 
"order by disorder"mL9 in which ground states that per- 



mit especially soft fluctuations about themselves are se- 
lected entropically at finite temperatures. 

This paper is concerned with the fate of such frustrated 
models, especially the maximally frustrated ones, when 
they are endowed with a quantum dynamics at T = 0. 
The canonical problems of this kind are S = 1/2 antifec-, 
romagnetic quanturp_Heisenberg models on triangularis 
and kagome latticesEj which have both attracted a con- 
siderable amount of interest. While one can think of these 
as quantized versions of classical Heisenberg models, it. 
is instructive instead to follow Anderson and FazekastS 
and think of them as perturbed variants of their Ising 
limits, which are manifestly frustrated. In these particu- 
lar cases, the perturbation is the XY exchange and as it 
does not commute with the Ising pieces, it introduces a 
quantum (two-spin) dynamics into the frustrated prob- 
lem. This perspective, in turn, suggests consideration of 
a more general class of models with alternative pertur- 
bations introducing a quantum dynamics instead. 

In this paper we report on a systematic study of frus- 
trated Ising models perturbed by transverse fields - as 
these introduce a single spin dynamics, thej4 are evidently 
the simplest models in the extended class.E-3 A shock par- 
tial, account of this work has appeared previously^ This 
paper covers a broader range of issues (and lattices) and 
supplies many of the details left out in the short account. 
It also provides a compendium of results we have ob- 
tained in the course of our work, which we hope will 
be of some use to people undertaking further study in 
this field. It is perhaps worth noting that this simplicity 
does make the models more tractable, which is gener- 
ally true for transverse field Ising models and accounts 
for their ubiquity in quantum statistical mechanical con- 
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texts. Much of this work is reviewed in Ref. |T^, which also 
includes some previous work on one-dimensional frus- 
trated chains. More recent studies include the use of 
the unfrustrated model as a paradigm of quantum crit- 
ical behaviorta and the treatment of random versions 
by asymptotically exact real space renonpalization tech- 
niques near infinite disorder fixed points J13 

A second motivation for studying these models is the 
possibility of direct experimental realization. Ising sys- 
tems exist and in the case of LiHoF4 and its kin, trans- 
verse fields have been used to, tune between phases in 
clean and disordered systems £3 A second family of Ising 
systems consists of stacked triangular lattices, strongly 
coupled along the stack, reviewed at length in Ref. ]LT\ . 
While these have an intimate connection to the single 
triangular lattice in a transverse field via its euclidean 
representation, it would be interesting in future work to 
consider the effect of a transverse field on them. Mean- 
while we would encourage our experimental colleagues to 
search for a triangular or kagome Ising antiferromagnet 
which are both, as we shall show in this paper, exem- 
plars of very different physics that can arise in frustrated 
quantum systems. 

Finally, there is the possibility of finding systems where 
there is a local Ising degree of freedom that is not itself an 
Ising spin. One such connection, which is currently the 
subject of intense interest, is an exacJ-jnapping from frus- 
trated transverse field Ising modelsE^I to quantum dimcr 
models of the short-range resonating valence bond state 
conjectured a while back by AndersonEj and oft men- 
tioned in discussions of the cuprate superconductors; in 
this context the Ising models aBpea|r more naturally as 
their dual, Ising gauge theories. E3~c3 

Turning now to the physics of the models themselves, 
we note that the introduction of quantum fluctuations 
can be expected to lead to a variety of behaviors much 
as in the case of thermal fluctutations catalogued above. 
Indeed, in the quantum case, the singular character of 
the fluctuations is manifest in that even infinitesimal 
strength perturbations lead to a non-trivial problem of 
degenerate perturbation theory in a macroscopically de- 
generate manifold and can therefore be expected to lift 
the degeneracy and select a much smaller number of 
ground states. Consequently there must be a discontinu- 
ity in the entropy and ground state correlations at zero 
quantum coupling and T = 0. We should emphasize that 
degenerate perturbation theory problems are not them- 
selves perturbative, especially for macroscopic degenera- 
cies. The most notable example of this is the quantum 
Hall problem where the degeneracy of a partially filled 
Landau level is lifted by the interaction and disorder in 
ways that lead to an incredibly complex phase diagram. 
Part of the interest of studying a diversity of perturba- 
tions of frustrated magnets is the prospect of generating 
at least a fraction of this complexity. 

Two possibilities are generic at small quantum rCOu- 
plings: a) a quantum version of order by disorderfl'EEl in 
which a broken symmetry state constructed out the the 



degenerate manifold is selected and b) a quantum version 
of the disordered or (cooperative) paramagnetic possibil- 
ity in which the ground state correlations remain short 
ranged and might be expected to lead to the opening 
of a gap in the spectrum. This latter possibility, which 
we have christened "disorder by disorder" ^ wag fast SU g_ 
gested by Anderson and Fazekascij as a mechanism for 
obtaining spin-liquid states of the RVB kind. Further 
singularities can emerge at large quantum couplings. For 
single spin dynamics, such as the transverse field prob- 
lems we study, there is neccessarily a paramagnetic phase 
at large couplings and at least one phase transition en 
route in cases of order by disorder. The nature of such 
phase transitions, which will be of some interest to us in 
this paper, can indeed be very unconventional. 
Specifically, we study Hamiltonians of the form: 



H — Jij Sf Sj 



(l.i) 



on a variety of one and two dimensional lattices (see 
Fig. ||) . Here the Jij are nearest- neighbor exchange cou- 
plings With \J zj \ = J and n p laqucttc(- J u/ J ) = ~^ T is 

the strength of the transverse field, the S a are the Pauli 
spin operators and h is the strength of a (classical) lon- 
gitudinal, Ising symmetry breaking field. 

The structure of this paper is as follows. We first in- 
troduce concepts and methods we have found useful in 
studying the model under consideration (Sects. || and 
HI). These we apply in the remainder of the paper to 
a number of frustrated Ising models on different lattices 
which, between them, realise a wide range of classical and 
quantum properties. We conclude by discussing possible 
realisations of such models. 

In detail, in Sect. ||, we use a mapping of the d dimen- 
sional quantum model to a d + 1 dimensional classical 
model to derive a general criterion for the existence of 
a quantum ordering transition, and we use this mapping 
to obtain a quantum Monte Carlo algorithm. This al- 
gorithm is free of any sign problems and is one of the 
attractions of studying this class of models - augmented 
by a cluster method appropriate for our problems, we 
have used i t at several points in this work. 

In Sect. Ill, we consider a three- leg ladder that has 
considerable pedagogical value and allows us to intro- 
duce a weak coupling (T/J <C 1) variational analysis 
and a strong coupling Landau- Ginzburg- Wilson analysis 
that will be our principal analytic tools in the remain- 
der of the paper. Sect. ^ reports results on the some- 
what baroque pentagonal lattice antiferromagnet which 
turns out to have an intimate connection with the lad- 
der described previously. Its classically disordered state 
gives way to quantum order. In the next three sections 
(Sect. [v|, [vi|, VII) we discuss three models: the trian- 
gular lattice antiferromagnet, the fully frustrated square 
lattice and the "two dimensional pyrochlore" lattice (see 
Fig. ||), respectively. These exhibit critical classical cor- 
relations in their ground state manifolds, and undergo 
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ordering transitions, in accordance with the ordering cri- 
terion derived in Sect. [n|. The latter also exhibits two 
unconventional critical phases. 

Next we turn to classically strongly disordered sys- 
tems. Sect. 



VIII 



deals with the kagome lattice antiferro- 
magnet which is a notable instance of disorder by disor- 
der and also exhibits a highly non-trivial phase diagram 
when both transverse and longitudinal fields are present. 
Sect. IX deals with a one dimensional Ising quantum dis- 
ordered magnet, the sawtooth chain, the classical ver- 
sion of which is the ultimate cooperative paramagnet. 
Finally, Sect. |x| analyzes the fully frustrated honeycomb 
lattice which appears to exhibit a fairly complex pattern 
of ordering driven by quantum fluctuations as well as a 
non-trivial, 0(4) phase transition. 

At various places in the paper we discuss connections 
to quantum height or dimer models. Via the latter, we 
find a connection to frustrated, valence bond, phases of 
Heisenberg magnets that .has been of interest starting 
from the oppposite end_J__ and-more recently starting 
from d-wave superconductorso'Ej As noted before, we 
are hopeful that this represents a more general possibil- 
ity of realizing the models we study here in other contexts 
where frustration is present and it is possible to focus on 
a local Ising degree of freedom. We close with a brief 
recapitulation of our themes in the summary section. 



II. A CRITERION FOR ORDERING AND A 
METHOD FOR QUANTUM MONTE CARLO 

In this section we map our e?-dimensional quan- 
tum magnet onto a d + 1-dimensional ferromagnetically 
stacked, classical magnet. This mapping allows use to 
derive some qualitative features of the quantum ordering 
process, and also to make contact with the existing lit- 
erature on stacked magnets. In addition, it enables us 
to devise a quantum Monte Carlo code free of the sign 
problem for general classical terms in the Hamiltonian. 



A. Mapping onto the stacked magnet 

We proceed j_.j_.by us- 

ing the Suzuki- Trotter formalism,c_'_3 to determine the 
precise correspondence between the quantum d and clas- 
sical d+1 dimensional models. Specifically, the partition 
function of the transverse-field Ising model is 



We now use a path integral representation of (2.1) 
where the insertion of a complete set of states effectively 
introduces an addi-t ional dimension of size /3, we follow 
Suzuki's approacI__: 
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Here, we have introduced the imaginary time step a T 
and n labels the coordinate of the extra dimension. The 
function 8^ k > is defined to be one if its arguments, the two 
(ordered) sets of spin configurations differ by k entries, 
and zero otherwise. 

We now establish an equivalence between the trans- 
verse field model and the classical Hamiltonian 

H d+1 = E K vV({Si, n }) + E K r Si, n Si, n+ i (2.4) 

(ij),n hn 

by expressing the partition sum for the latter in terms of 
a transfer matrix, Tz'. 



T z = exp(K v V({S itn })) x 
+ 0(exp(-iT))) . 



(2.5) 



The first term on the right hand side of the previous 
equation is to be understood as an exponentiated diago- 
nal matrix. 

^From this, one can read off that the two partition 
functions will be equivalent if one chooses a r r/2 = 



exp(-K T /2) and K{ 



We note that continu- 



ous quantum evolution corresponds to the scaling limit 
Ky oc a T — » 0, K T — » oo, while maintaining 



2e 



/a r = r, K v exp(K T /2)/2 = v/T 



(2.6) 



Z = Tre-i 3H = Trexp <^ -0 W{{Sf}) 




(2 



where = 1/ksT. For generality, we have introduced the 
notation V({S*}) for the 'classical' part of the Hamilto- 
nian of strength v, i.e., [V, Sf ] — for all i; in the sim- 
plest case, V only consists of the exchange part of H, 



With the classical Hamiltonian H^+i (Eq. ^4|), the 
Ising spins interact in the spatial layers as they do in the 
analogous d dimensional classical problem, but they are 
l^also coupled ferromagnetically in the additional (imagi- 
nary time) dimension. The dimensionless inverse of the 
quantum temperature, 0T, is given by the extent L T 
of the system in the time direction: 0T — a T TL T = 
2exp(-K T /2)L T . 
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B. A criterion for order by disorder 

We can look at the possibility of quantum ordering in 
a transverse field by studying the discretized partition 
function Z OQ (K T ) = Z(K S — > oo, K T ) as an expansion in 
powers of K T . In this limit, we force the spin configura- 
tions in each plane to be classical ground states. Taking 
the trace, Tr gs over these ground states gives 



Zoo(K T ) =Tr gs exp K T 2J Si, n Si, n +i 



(2.7) 



where the sum is over all ground-states for each layer. 
Expanding in small K T , we obtain 



1 



Z oc {K T )= Tr gs ( 1 + 1(K T ) 2 x 



x /J y ^ i.n^-t.n+l^n'^j.n'+l + 0((K T ) 4 

where the linear term in K T is abse nt d ue to Ising sym- 

as 



metry. We can further express Eq. (2 
Z X (K T ) = Z{L) L 



|i + i(^ r ) 2 E<^„^„ +1 ) 2 + ...| 



where the prefactor refers to the number of classical 
states in each layer of area 1? and L T = 0/a T is the 
number of layers. Then 



Zoo(K T ) = exp(L 2 L T 5) x 



(2.10) 



|l + i(A-) 2 (L 2 ^)E(^ 5o,o) 



where S is the classical ground-state entropy density. 
Therefore the effective free energy, i 7 ^ = — lniv 00 //3 as 
a function of small K T is 



(3F X = (L 2 L T ) { S 



{K 7 



(2.11) 



where the term (L 2 L T ) is a volume in space-time. 

This yields the following powerful result. If the sum 
/ = J2i(SiflSo,o) 2 °f the classical correlation function di- 
verges, the free energy above is non-analytic as K T — > 0, 
implying that the K?j = oo is in a different phase than 
the disordered point K T = for any K T > 0; this is 
indeed the case for the Ising triangular antiferromag-j 
net where Itri contains a leading divergence of the formcl 
/ d 2 r(l/v^) . In this case, one does indeed find quan- 
tum order by disorder (see Sect. 0). 

However the situation is inconclusive for classically 
disordered antiferromagnets with exponentially decaying 



spin correlations at T = ; eve n though no single term in 
the series expansion, Eq. ( |2.11 ) , becomes unbounded, the 
whole series may diverge because it is beyond its radius 
of convergence. There is therefore no clear link between 
classical and quantum disorder. 



C. Quantum Monte Carlo 

By simulating the stacked classical magnet, it is thus 
possible to gain information on the properties of the 
quantum system. We note that, thanks to the simplicity 
of our quantum dynamics, there is no sign problem to 
cope with for this class of models. 

However, the scaling limit K T — ► oo does pose some 
technical problems. At large K T , domain walls in the 
time direction become very rare, which leads to a di- 
vergent timescale in the Monte Carlo simulations. This 
problem can be remedied by employing a cluster algo- 
rithm, in which the attempted Monte Carlo moves con- 
sist not of flipping a single spin but rather rods of spins 
in the time direction. These rods can be chosen in a 
way which exactUz-.cancels the inclement Boltzmann fac- 
tor exp(— K T /2)E3 It is with this method that we have 
carried out the simulations presented here, with an ad- 
ditional feature for the fully frustrated hexagonal lattice, 
(^which is described there (see Sect. g). 

It however turns out that the dominant source of error 
can be the remaining discretisation error. To see this, 
consider the case of a quantum ordered state at zero 
quantum temperature, 1/(3 = 0, which corresponds to 
an infinite extension of the stacked model in the time di- 
rection. One can now imagine starting at an effectively 
infinite classical coupling, K T . As K T is reduced, within 
the framework of the above mapping, one retains 1//3 = 
but a discretisation error is introduced - which for suf- 
ficiently small K T will be serious enough to make the 
quantum order disappear as the classical magnet goes 
through its transition at a critical K T . If one consid- 
ers a finite /3, one in practice has to trade off this error, 
which requires large K T against reaching a low quantum 
temperature, which nominally decreases with increasing 
K T . 

To quantify the discretisation error, we quote the pa- 
rameter A = exp(kK T /2), (where fc is the multiplicity 
of the spin flip, being one for single-spin dynamics), a 
lengthscale characteristic of an isolated ferromagnetically 
coupled rod - the larger A, keeping X/L T fixed, the better. 
For fixed L T , however, the optimal value of A depends on 
the correlation length in the time direction, which varies 
from system to system and which is larger the smaller 
the quantum gap. 

In addition, the time for building up the clusters also 
becomes large at low temperature. This reflects the fact 
that the stacked classical magnet represents a rather inef- 
ficient way of doing the bookkeeping for the spin state. In 
the presence of only a few domain walls, it is probably su- 
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perior to keep track of the domain walls themselves. This 
can be done in the framework of a continuous time algo- 
rithm, recommended for future use, which is described 
for example in Ref. [52]. 

III. METHODS USED AND THE FULLY 
FRUSTRATED THREE-LEG LADDER 

To introduce some of the concepts used repeatedly in 
this paper, let us first consider the toy model of the fully 
frustrated three-leg ladder depicted in Fig. ||. The inter- 
actions along the rungs are ferromagnetic, as are those 
along the outer two legs. The antiferromagnetic inter- 
actions along the middle leg make the ladder fully frus- 
trated. 

The classical ground states are those states which min- 
imise the number of frustrated bonds. Since the bonds 
on the outer two legs belong to only one plaquette, it 
is not favourable to frustrate these. By contrast, frus- 
trating one inner bond can put the two plaquettes it be- 
longs to intOj-the ground state. One finds that there are 
three sectorscJ of ground states, and their Ising reversed 
counterparts, depicted in Fig. g. In all these states, the 
top and bottom legs are ordered ferromagnetically. The 
staggered sectors, with the spins on the top and bottom 
legs of opposite sign, contain only one ground state each; 
by contrast, the columnar sector has an extensive zero- 
temperature entropy per rung of S = So/^Nks) = G, 
where G = {\fE + l)/2 is the golden mean, which can be 
obtained by a transfer matrix approach. 

Since each plaquette has exactly one frustrated bond, 
which is shared with a neighbouring plaquette, one can 
represent the ground states by placing dimers onto the 
frustrated bonds. This gives rise to a hardcore dimer 
covering of the dual lattice, which in this case is a two- 
leg ladder. Note that such a mapping is possible for all 
fully frustrated models in which the elementary plaque- 
ttes are arranged to share bonds. It does not, however, 
exist for all lattices; for instance, if the plaquettes share 
sites rather than bonds, as is the case for the kagome lat- 
tice, the hardcore nature of the dimer model is lost. The 
advantage of this mapping is that, up to a global spin re- 
versal, there is a one-to-one correspondence between the 
ground states and the hardcore dimer coverings. There- 
fore, restricting the full Hilbcrt space to that of the dimer 
coverings yields a natural implementation the projection 
onto the ground state. Below, we will present a deriva- 
tion of the transverse field Hamiltonian restricted to the 
dimer (ground-state) manifold. 



A. The action of the transverse held — mapping onto 
a quantum dimer Hamiltonian 

In the S z basis, the transverse field operator TS X is 

a spin-flip operator, (r/2) ^ ^J. The leading-order 

quantum dynamics within the ground state manifold is 
therefore that of single spin flips connecting different 
ground states - we need to consider multiple spin flips 
only in cases where single spin flips are nowhere possi- 
ble, as will be the case in some of the examples discussed 
further down. 

Flipping a spin without leaving the ground state mani- 
fold is possible only if the spin is part of the same number 
of satisfied bonds as frustrated ones, or, in other words, 
if it experiences zero net exchange field. We refer to such 
spins as flippable spins. In dimer language, such a spin 
is at the centre of a dimer plaquette of the form I I or — , 
and flipping the spin exchanges the frustrated and the 
satisfied bonds and thus corresponds to the elementary 
dimer move I I <-> — . 

The two staggered sectors of the ground state contain 
only the two staggered dimer states. The states in the 
columnar sector can be obtained starting from the colum- 
nar dimer state by repeated application of the elementary 
plaquette move I 1 — ► ~. 

Therefore, in the staggered ground-state sectors the 
transverse field has no effect since no elementary dimer 
moves are possible - indeed, there are no degenerate 
states which it can mix at T = at any finite order 
of the perturbation theory. This is different in the main 
sector, where the transverse field lifts the macroscopic de- 
generacy and promotes a particular linear combination of 
the classical ground states to the true, quantum ground 
state. 

In fact, it is now apparent that the transverse field 
Hamiltonian, restricted to the classical ground state man- 
ifold, in dimer language can be written as 

Hqdm = ~t (|~> (I I| + h.c) + v (p) p + |II)(I I|) , (3.1) 

where the kinetic term with t = r/2 is generated by the 
transverse field. We have added a diagonal term with 
coefficient i>, which is zero for the transverse field problem 
but which will be useful later on. 

Note that this Hamiltonian has only nonpositive off- 
diagonal matrix elements, so that the Perron-Frobenius 
theorem can be used to predict a nodeless quantum 
ground state, which means that in it, the amplitudes 
of all the configurations can be chosen to be real and 
nonnegative. 

Since a ground state and its Ising reversed counter- 
part both get mapped onto the same dimer states, we 
have to show that the dimer ground state, is in fact the 
same as the transverse field ground stated First we split 
the ground state manifold into two submanifolds, namely 
those containing symmetric and antisymmetric combina- 
tions of Ising-reversed pairs of states. The transverse field 
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Ising Hamiltonian does not connect these submanifolds, 
so that it is block-diagonal. 

All entries in the symmetric block continue to be of 
the same (negative) sign so that the ground state in this 
block continues to be nodeless. Since the full Hamilto- 
nian is block-diagonal, the state obtained by combining 
this ground state with a null state in the antisymmet- 
ric block continues to be a nodeless eigenstate of the full 
problem. Since there is only one nodeless eigenstate, the 
state obtained by translating the dimer ground state into 
the spin ground state is the ground state of the transverse 
field problem. 

The entries in the antisymmetric sector have the same 
modulus as in the symmetric one; however, they need 
not all have the same sign. If we pick one spin state 
for each dimer state and collect those in the up mani- 
fold, and their spin-reversed counterparts into the down 
manifold, an entry will be negative when the transverse 
field connects members of the up and down manifolds. 
The states can be sorted in a way, e.g. by magnetisa- 
tion, that the fraction of negative entries vanishes in the 
thermodynamic limit. Whether the ground states in the 
two sectors in the thermodynamic limit are degenerate 
then depends on whether the wavefunction has substan- 
tial support on the states at zero magnetisation, in which 
case they are not, or whether it is localised away from 
them, in which case they are. The example of the three- 
leg ladder is special in that the up and down sectors can 
be chosen to be entirely disconnected, so that an exact 
degeneracy trivially arises. 

The connection between transverse field Ising model 
and quantum dimer model we have established is use- 
ful for several reasons. It affords some insight into the 
structure of the problem we are studying in that it pro- 
vides a natural deformation of the transverse field Ising 
model by switching on the potential term, i.e. by choos- 
ing a nonzero v. The cases \v\ ^> t are easily solved and 
can therefore provide two anchors of the phase diagram 
containing the point we are interested in. 

Moreover, there is the special point v — t, known 
as Rokshar-Kivelson (RK) point, after the inventors 
of the modelpa where the (nodeless) quantum ground 
state is an equal-amplitude superposition of all classi- 
cal ground states. Therefore, operators which are diago- 
nal in the dimer basis provide precisely the expectation 
values of the corresponding classical operators. In this 
spirit, switching on an infinitesimal transverse field im- 
plies jumping the finite distance from v = t to v = 
in the Rokshar-Kivelson model (Fig. ||), a manifestation 
of the fact that the transverse field perturbation is non- 
analytic. Deciding the ordering pattern of the transverse 
field problem can therefore be accomplished if one can 
show that the RK and the transverse field points are in 
the same phase. 

Moreover, the transverse field problem in itself pro- 
vides a new perspective on the Rokshar-Kivelson quan- 
tum dimer model, which was proposed as a model of 
Anderson's (short-range) resonating valence bond (RVB) 



physicsBIl It can be derived for a Heisenberg antiferro- 
magnet: the perturbative derivation uses the nonzero 
overlap between different dimer configurations as expan- 
sion parameter. This model is useful in a regime where 
the Heisenberg model is in a phase dominated by valence 
bonds. The study of a transverse field Ising model on a 
frustrated lattice can thus be used to gain insight into 
the behaviour of Heisenberg magnets on the dual lat- 
tice. This fact has been used by the present authors to 
identify a bonaJide short-range RVB phase on the tri- 
angular lattice,E3 a goal that had proven to be elusive 
on the square lattice, for which the model was originally 
formulated. 

Finally, we note parenthetically that for the mixed 
bond models, the application of a longitudinal field is 
somewhat arbitrary as there is a (gauge/Mattis) freedom 
of which bonds to call antiferromagnetic and which ferro- 
magnetic, as long as the odd condition is met. Thus any 
state in the gauge invariant dimer representation can be 
(up to topological restrictions) be represented by a max- 
imally polarised (long-range ordered) spin configuration, 
irrespective of the nature of the dimer correlations. In 
the models where all bonds can be chosen to be anti- 
ferromagnetic, this choice defines the natural gauge and 
thus makes the application of a longitudinal magnetic 
field unambiguous. 

B. Maximally flippable states 

We now present a heuristic argument which is varia- 
tional in spirit to generate a candidate state for selection 
by the transverse field. This state is the one which can 
gain the most energy from the transverse field on account 
of being composed around a backbone configuration hav- 
ing the strongest fluctuations possible. 

To start, we note that in order to gain energy from the 
transverse field, spins have to have a component pointing 
in the ^-direction, which in the S^-basis goes along with 
the component of the form [| T) + I I)] /V%- Those states 
in which most of the spins can be flipped and hence put in 
a superposition of up and down states thus stand to gain 
most energy from the fluctuations induced by the trans- 
verse field. We therefore identify the configuration with 
the largest number of flippable spins ("maximally flip- 
pable state") as the backbone of our candidate ground 
state. The actual ground state naturally includes fluc- 
tuations around the maximally flippable configuration, 
since it,is these which lead to the energy gain in the first 
place.Ej The energy gain due to alignment of the spins 
along the transverse field in the quantum model is an 
entropic contribution to the free energy in the classical, 
stacked model. 

As an aside, we note that the classical correlations pro- 
vide the simplest first guess at the maximally flippable 
configurations and hence at the ordering pattern. This is 
because configurations with many flippable spins have a 
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large number of neighbouring configurations which differ 
only by a few spin flips. The correlations they incor- 
porate, even when not leading to long-range order, can 
thus already be visible in the classical average where all 
ground states are accorded equal weight. 

For our three-leg ladder, the maximally flippable con- 
figuration is the columnar one depicted in Fig. ||, since 
there, all spins on the middle row (all of which point up) 
are flippable. We therefore expect to find a state which 
has a ferromagnetic moment even in the z-direction in 
addition to the polarisation in the re-direction. 

The selected state can incorporate structure in addi- 
tion to that apparent from the maximally flippable con- 
figuration. This follows from the fact that although the 
backbone configuration maximises the number of flip- 
pable spins, in fact not all the spins are independently 
flippable. In the three-leg ladder, for example, flipping a 
spin on the middle leg precludes flipping its neighbours, 
so that in effect only half the spins are independently 
flippable. One can therefore give two different recipes for 
constructing the quantum state. 

To contruct the first type of state, we take all the flip- 
pable spins and polarise them in the x-direction disre- 
garding the ground-state constraint. Next, we reinstate 
the ground-state condition by projecting out those com- 
ponents of the state which are not contained in the classi- 
cal ground-state manifold. We call this state the uniform 
state since it treats all the flippable spins on the same 
footing. 

The second type of state is obtained by identifying the 
largest set of independently flippable spins and polarising 
those in the rr-direction. Since, starting from a maximally 
flippable configuration, there are typically several choices 
for which set of flippable spins to polarise, we call this 
state the hierarchical state. It will have a lower symme- 
try than the maximally flippable configuration unless all 
flippable spins are independently flippable. 

The distinction between the uniform and the hierar- 
chical states will turn up in several times in this article, 
and it arises quite naturally in other approaches. Note 
that this approach suggests yet another class of alter- 
native candidate configurations, namely those which do 
not have the maximum number of flippable spins but 
nonetheless maximise the number of spins which are flip- 
pable independently. 

The variational states thus obtained have the short- 
coming that all their component configurations are allo- 
cated equal weight. This is clearly not optimal, since the 
ground state will in any case place maximal weight on the 
individual, maximally flippable configurations. Within a 
more elaborate variational framework, weights could be 
accorded to the configurations depending on their num- 
ber of flippable spins. 

In summary, the flippability approach suggested here 
identifies a 'saddle point' - the maximally flippable con- 
figuration - that will be favoured due to the fluctuations 
around it, which are 'softer' than those around other con- 
figurations. We expect it to work as long as the actual 



quantum wavefunction is concentrated on the maximally 
flippable and nearby configurations. It will break down 
if the wavefunction has the bulk of its support elsewhere, 
i.e. on configurations unrelated to the maximally flip- 
pable one. 

This point can be made more intuitively by considering 
the degenerate perturbation theory as a hopping prob- 
lem. Each classical ground state configuration defines a 
point in the ground-state manifold. The transverse field, 
by flipping spins, connects different configurations, thus 
endowing the ground-state manifold with a graph struc- 
ture. The perturbation theory can be thought of as a 
hopping problem on the graph thus defined. The large 
weight of maximally flippable configurations follows from 
their high coordination, and the magnetic ordering tran- 
sition corresponds to a localisation transition in the hop- 
ping problem E3 

These ways of thinking are in close correspondence to 
the case of thermal order by disorder n There, thermal 
fluctuations (out of the ground state manifold) provide a 
large entropic weighting to the states allowing the softest 
fluctuations. These states are then selected as T — > + 
provided their enhanced weight is not swamped by the 
combined fluctuatioHal and configurational entropy of the 
less favoured statesO However, when the thermal fluctu- 
ations increase in strength with increasing temperature, 
they destroy the order which they were instrumental in 
establishing in the first place, as happens for T large in 
our problem, as discussed in the following paragraphs. 

C. The opposite limit: Landau-Ginzburg analysis for 

r> j 

In addition to doing the (hard) degenerate perturba- 
tion theory for small transverse fields, we can use an alter- 
native approach for determining the state of the quantum 
magnet which is made tractable by virtue of the simplic- 
ity of the transverse-field term. Consider the problem 
where the relative sizes of exchange and transverse field 
are inverted, namely where r> J. In the limit J/T = 0, 
the ground state is a simple paramagnet in the S z basis: 
all spins are polarised along the -(-^-direction explicitly 
selected by the field. In addition, it is gapped: the low- 
est excitations are spin flips which each cost an energy of 
T; this makes it possible to perturb about this state by 
switching on a weak exchange. This contrasts to the case 
where the transverse field is replaced by an XY exchange. 
Here, the large Jxy problem is not exactly soluble, and 
so it is not possible to perturb about it - indeed, it may 
even be gapless. 

As the perturbing J is switched on in addition to the 
transverse field, the excitations acquire a dispersion, typi- 
cally - but not always - already to first order in J/T. The 
dispersion to first order is simply given by the Fourier 
transform of the interaction matrix of the lattice. For 
sufficiently large J, there can be an ordering transition, 
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which corresponds to a macroscopic occupation of the 
softest mode(s). 

To generate the state to which the leading transition 
takes place, one has to combine knowledge of the soft 
modes with lattice symmetry considerations to construct 
a Landau-Ginzburg- Wilson action. This program is car- 
ried out in detail for the fully-frustrated hexagonal mag- 
net below and follows the wack on layered magnets by 
Blankschtein and coworkers. l^Ilj 

The combinations of the soft modes dictated by sym- 
metry considerations yield the ordering pattern which 
is established as J/T is increased. In addition, by also 
analysing the Landau-Ginzburg- Wilson action with stan- 
dard methods of the renormalisation group, one can ob- 
tain information on the nature of the transition into the 
ordered configuration, and it turns out that these tran- 
sitions are generally not Ising transitions, as the most 
naive guess would suggest. By analysing the potential 
presence of dangerously irrelevant terms in the action, 
one can even guess at further symmetry breaking. 

This analysis of course has the usual limits associated 
with mean-field theories. Most important in this context 
is the possibility of further phase transitions out of the 
ordered phase before we reach the regime of infinitesimal 
transverse fields, where the 'small' parameter J/T — > oo. 
Another scenario is the absence of any phase transition, 
so that the magnet remains disordered at all couplings. 
Although the excitation dispersion may soften at par- 
ticular points in the Brioullin zone, fluctuations may be 
sufficiently strong to prevent ordering at any coupling. 
However, even in this situation, the large T approach 
may be used successfully if the expansion in powers of 
J/T is carried to sufficiently high order and combined 
with a non-perturbative analysis such as that given by 
the use of Pade approximants. In Ref. |40| this program 
has been carried out for the sawtooth chain (see Fig. ||). 

It is worth pointing out that this approach presents, 
ultimately, a soft-spin analysis in that the size of the or- 
dered moment can vary from site to site so that, as the 
mode amplitude increases, non-linearities become impor- 
tant. However, even for our hard Ising spins, a difference 
in the size of the ordered moment does have a meaning. 
As an illustration, consider the uniform state of the three- 
leg ladder, defined above. Here, every spin on the middle 
row is fluctuating and thus has \(S Z )\ < 1/2, whereas 
the remaining spins have (S z )=l/2, being fully polarised 
along the z-axis. This difference is indeed found in the 
large T analysis and can thus be interpreted as being due 
to the fluctuations induced by the transverse field. 

In detail, for the three-leg ladder, the Fourier trans- 
form of the interaction matrix is (labeling the sites on 
the rung 1 to 3 from top to bottom) 

T / 2 cos(fc) 1 \ 

- 1 -2cos(fc) 1 , (3.2) 

1 \ 1 2 cos(fc) / 

from which we obtain a degeneracy in the dispersion re- 
lation. One minimal mode, with wavevector q = 0, has 



eigenvector (1,-2 + v6j 1); the other has q = it and 
eigenvector (1,-2 — y/E, 1). Here, q is the wavevector 
along the ladder, and the entries in the eigenvector de- 
note the amplitudes of the top, middle and bottom sites 
(\/6 ~ 2.45). 

This degeneracy is accidental in that the states are not 
related by any symmetry operations, and we have not en- 
countered this effect in any of the more regular lattices 
we are studying. Quite generally, however, in lattices 
with sites of different coordination, the soft-mode analy- 
sis might suggest states which fare very poorly under the 
hard-spin constraint, as is the case for the q — ir mode 
in this example (see below). A more appropriate order- 
ing pattern may nonetheless show up as an alternative, 
possibly only local minimum in the dispersion relation. 

The q = mode can easily be identified as the actual 
ordering pattern (see below, exact diagonalisation) , with 
a moment on the middle row reduced by the fluctuations. 
One reason the competing state loses out eventually is be- 
cause it has a reduced moment on the sites which are in 
fact not allowed to fluctuate in the small-r limit: there, 
the leading term in the perturbation theory flips spins 
with an equal number of frustrated and satisfied bonds, 
and the reduced moment sites in this state have odd co- 
ordination and are thus never flippable. 

D. Exact diagonalisation 

For the three-leg ladder, one can attack the transverse 
field problem by exact diagonalisation of the degener- 
ate perturbation theory. Systems with up to fourteen 
rungs are easily accessible numerically. The results are 
depicted in Fig. |[ We do indeed find a state with ferro- 
magnetic order along the centre row, as predicted by the 
flippability analysis. The selected state is the uniform 
columnar one, as the gap extrapolates to a finite value 
for large system sizes and hence the breaking of transla- 
tional symmetry accompanying the hierarchical state is 
absent. 

The infinitesimal transverse field generates a nonzero 
polarisation in the x-direction, corresponding to an 
ground-state energy per spin of E = — (0.302 ± 0.001) |T| 
in excess of the classical value. The finite moment (S z ) al- 
ready present in the classical ground state average (where 
it equals 1/(2(2G— 1)) ~ 0.224), is discontinuously en- 
hanced as r is switched on. As the system size is in- 
creased, the maximally flippable state and its neighbours 
rapidly gain in weight, whereas at the other extreme, 
those states without flippable spins have zero amplitude 
in the ground state wavefunction. 

These results tie in with the picture provided by the 
mapping onto a quantum dimer model. The selected 
columnar dimer phase, which extends all the way to 
v = — oo, terminates at the RK point, beyond which 
it gives-way to the staggered dimer phase, as depicted in 

Fig. pa 
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E. Relation to XY perturbations 

Whereas the transverse field induces single spin flips, 
an XY exchange flips ncigbouring pairs of antialigncd 
spins. For lattices with even coordination and all bonds 
antiferromagnetic, the XY flips are a subset of the two 
spin flips generated by two successive operations of the 
transverse field. 

For lattices with odd coordination, however, the trans- 
verse field can never flip a spin without leaving the 
ground state manifold, so that the XY moves can show 
up as the leading terms, in order T 2 / J, in the degenerate 
perturbation theory. These matrix elements are negative, 
corresponding to a ferromagnetic XY exchange. This will 
be of importance for the fully frustrated hexagonal lat- 
tice. Moreover, in the case of the pentagonal lattice, 
where inequivalent sites with even and odd coordination 
exist, the moves generated by an XY perturbation are 
entirely distinct from those generated by the transverse 
field. 



IV. THE PENTAGONAL LATTICE 

The pentagonal lattice, depicted in Fig. ||, tiles the 
plane with (irregular) pentagons. It can be obtained from 
the hexagonal lattice by cutting each hexagon in half with 
a set of parallel lines ('cuts'). Classically, the Ising anti- 
ferromagnet on this lattice is disordered and has ariinite 
ground-state entropy per spin of S/ks = 0.234.113 The 
ground states can be represented by a hardcore dimer 
model as described above. 

The flippability analysis for this magnet is straight- 
forward. To lowest order in T, only the spins along the 
cuts can be flipped since only they have even coordina- 
tion. Therefore, the spins along the sawtooths follow an 
antiferromagnetic pattern, since the maximally flippable 
configuration (depicted in Fig. ^ has all frustrated bonds 
associated with spins on the cuts. This configuration is 
closely related to the fully frustrated three-leg ladder in 
that the spins along the cuts are effectively decoupled 
from the other spins in the system and the transverse 
field generates the same Hamiltonian as for the ladder. 
Thence, the ordering pattern along the cuts will be the 
q = pattern described above, with a q — ir modulation 
transverse to the cuts. 

For completeness, we now perform the soft-mode anal- 
ysis for this problem. The interaction matrix is 





There is a line where the eigenvalues attain the global 
minimum — J(l + VS)/2 at wavectors (tt, y). The eigen- 



vectors ((—I — y/S)e 
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, l) correspond to states 



to lowest order in the degenerate perturbation theory 
in r, and this state is therefore not competitive in the 
small r limit. There is, in addition, one further local 
minimum — Ja/6/2 at wavevector (0,7r). The eigenvec- 
tor (—2 + \/6, — 1, l) corresponds to the maximally flip- 
pable state (and agrees with the one given by the three- 
leg ladder). We therefore expect a phase diagram with 
at least three transitions coming from large T, namely 
first from the disordered into the 'sawtooth state' fol- 
lowed by a transition into the maximally flippable state, 
which finally terminates in the disordered classical phase 
at r = 0. There may be a separate transition which se- 
lects one or a combination of the sawtooth states - we 
have not investigated this. The complicated structure of 
this phase diagram is of course to a large degree a conse- 
quence of the presence of inequivalent sites. 



V. THE TRIANGULAR ANTIFERROMAGNET 

The two-dimensional Ising antiferromagnet on the tri- 
angular lattice has power-law spin correlations and an 
extensive entropy (S/k B = 0.323) at T = OBBtf We first 
describe a number of useful mappings and then discuss 
their implications. 



A. Mappings to height and dimer models 

In Fig. we have displayed a particular spin configura- 
tion. The corrcsonding dimer configuration (also shown) 
is constructed as in Sect. III). The effect of the trans- 



with frustrated bonds on the sawtooths. As mentioned 
above, these do not lead to flippable spin arrangements 



verse field can easily be stated in terms of the dimer map- 
ping in a more general setting than the one discussed in 
Sect. III). If we choose a closed loop of bonds on the 
dual lattice which are alternately empty and occupied 
by a dimer, and exchange the occupied with the empty 
bonds, this corresponds to flipping all the spins inside 
the closed loop. The leading order effect of the trans- 
verse field is therefore to produce such a rearrangement 
for the loop enclosing the smallest number of sites, and 
the relevant order in perturbation theory is given by the 
number of enclosed spins. The loop in question is typi- 
cally the shortest closed loop of even length on the dual 
lattice. In this case, this is a loop of length six enclosing 
one site (see Fig. [7]), but we will encounter the case of 
more than one enclosed site later on. 

Next, we turn our attention to the mapping to a height 
model. The height configuration is given by the numbers 
in Fig. 0. The height variables are defined on the sites of 
the lattice.EJ Height differences are determined as follows. 
If one crosses a (no) dimer when going along a bond 
in a clockwise direction around a triangle pointing up, 
the height changes by 2 (—1). If one goes the opposite 
direction, the height change is —2 (+1). One can easily 
check that this prescription has the property that the 
height change around the three sides of any triangle is 







zero. Moreover, any closed path can be decomposed into 
a sum of paths around individual triangles, and thus the 
height change between any two sites is indepent of the 
path chosen to go between them. 

The existence of the height model has important ram- 
ifications for the classical model. Under a set of reason- 
able assumptions).^ the result follows that the ground 
state correlations are algebraic. Combined with our cri- 
terion for quantum ordering (Sect. II B), one expects all 
magnets allowing a height mapping to order. The order- 
ing correlations are then determined by the configuration 
which is as flat (in height space) as possible. 

This raises the question under which conditions a 
height mapping exists. A set of sufficient conditions is 
that, firstly, the lattice be composed of bond sharing frus- 
trated units (allowing a dimer model) and secondly, its 
dual lattice be bipartite. This is required for giving a 
set of consistent rules for computing height differences. 
In the cases we encounter below where these conditions 
are not met, we will see that a height mapping is absent, 
consistent with their disordered classical correlations. 



B. Ordering behaviour of the triangular IAFM 

We first study this problem for small transverse field 
with h = 0. The maximally flippablc configuration is the 
one depicted in Fig. |?|. It has a unit cell of three sites 
(two spins pointing up and one down). This is also the 
flat (height) configuration, as well as being the state en- 
countere|d at v = — oo in the appropriate quantum dimer 
model:E3UJ 



H, 



QDM = -t (|v>(£| + h .c) + V (|V)(V| + 

In fact, this state has a net magnetisation since its 
three-sublattice strucutre is (1,1,-1), with all spins 
pointing up being flippable. The uniform state will thus 
have a sublattice structure of the form (a,a,—b), with 
b > a, whereas the hierarchical one will have (c, 0, — c). 

All the evidence therefore points towards a three 
sublattice quantum ordering pattern, with details to 
be determined. Such correlations were previously 
found by a Landau-Ginzburg- Wilson (LGW) r afialysis of 
ferromagnetically-stacked triangular lattices]^ applica- 
ble to our problem by virtue of the mapping presented in 
Sect. II A. This LGW analysis finds an action which is 
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XY symmetric up to sixth order, where an XY symmetry 
breaking, six-fold clock term appears. 

The resulting phase diagram in the temperature-field 
(T — r) plane, depicted in Fig. ||, is quite remarkable. At 
T = 0, the triangular TFIM undergoes a quantum phasa 
transition which is in the 3-d XY universality class,Ej 
where the clock term is dangerously irrelevant, so that 
the transition is into a phase with clock-symmetry break- 
ing as well. 

At finite temperature, however, the transverse field 
Ising model maps onto a classical stacked magnet of fi- 
nite size, L T , in the imaginary time direction, so that 



its behaviour crosses over to being two-dimensional as 
the correlation length becomes comparable to L T ; in this 
regime, we thus have to consider the properties of a two- 
dimejasional XY model with a clock term. Here one 
findsc3 two finite-temperature Kosterlitz-Thouless tran- 
sitions bordering an extended critical phase. The phase 
diagram in Fig. @ is reliable in the region near the three 
dimensional XY transition, where the order parameters 
are small. 

Depending on the sign of the six-fold clock term, the 
quantum ordered state is predicted to have the sublat- 
tice structure of either the uniform or the hierarchical 
state. In order to confirm the predicted ordering pat- 
tern and to find out the precise sublattice structure, we 
have carried out quantum Monte Carlo simulations on 
the triangular TFIM. The results are depicted in Figs. ^ 
and 10. A Bragg peak is clearly visible at q x — 2tt/3, as 
expected. The sublattice structure of the form (1, 0, —1) 
is also clearly borne out, showing that the hierarchical 
state without a net moment is selected. Since the quan- 
tum dimer state at v = — oo is the uniform, magnetised 
state, this implies a transition for v < 0. Preliminary 
studies by us indicate that such a transition is indeed 
present close to v = 0. Finally, we note that a finite 
longitudinal field (\h\ < 6J) on its own selects the same 
uniform state. 



VI. VILLAIN'S ODD MODEL: THE FULLY 
FRUSTRATED SQUARE LATTICE 

A great deal is known about the square FFIM, espe- 
cially when we combine the knowledge about the differ- 
ent models it is equivalent to by virtue of the mappings 
presented above. It turns out to shadow the triangular 
IAFM in practically all important qualitative respects. 
Its classical ground-state entropy S — C/n ~ 0.241,eZI 
where C is Catalan's constant. Since it consists of bond- 
sharing squares with a bipartite dual lattice, it is critical 
at T = 0. The cor responding dimer model Hamiltonian 
is given by Eq. 3.1 and has a phase diagram like that pic- 
tured in Fig. a. The columnar phase corresponds to the 
maximally flippable state which, in dimer language, is ob- 
tained by covering the square lattice with infinite dimer 
ladders with the configuration depicted on the right panel 
of Fig. [|. An LGW analysisEj again finds an XY ac- 
tion with an XY symmetry breaking clock anisotropy at 
higher order, which, depending on its sign, selects either 
the uniform or the hierarchical (also known as plaquette) 
state. 

This makes it clear that there is quantum ordering 
in the FFIM on the square lattice into a state with 
translational symmetry breaking. We have not stud- 
ied this problem numerically ourselves but we refer the 
reader tGirtte literature on the ferromagnetkrally stacked 
magnetEjO and the quantum dimer modeloo 
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VII. PYROCHLORE IN TWO DIMENSIONS: 
THE SQUARE LATTICE WITH CROSSINGS 

This lattice, also known as the checkerboard or two- 
dimensional pyrochlore lattice, is made up of a square 
lattice with nearest neighbour interactions which has, in 
addition, crossing next-nearest neighbour interactions on 
alternate plaquettes as depicted in Fig. [ll]. The plaque- 
ttes with crossings are equivalent to tetrahedra in that 
they contain four sites all interacting equally with one an- 
other. Since these tetrahedra are arranged to share sites, 
as is the case in the pyrochlore lattice in three dimensions, 
from which it can in fact be obtained by a projection in 
a (lll)-direction; for a picture, see Ref. |37|b. 

The ground-state condition for the classical Ising 
model is that each tetrahedron have zero magnetisation. 
There are six such states for each tetrahedron, each with 
two spins up and two down. The Ising model an the 
pyrochlore lattice is equivalent to the ice modelpJ the 
ground state entropy of which was calculated exactly 
in two dimensions by Lieb:Lj S/ks — § In | w 0.216. 
This model can be mapped onto a six- vertex model (each 
vertex encoding one of the single-tetrahedron ground 
states), and hence onto a non- intersecting loop model, 
which in two dimensions guarantees the existence of a 
height model. As discussed above, this quite gener- 
ally implies classically critical gijound-state correlations, 
which have indeed been found £3 

Here, we present a short derivation of the height model 
which we then use to determine the quantum ordering 
behaviour; refer to Fig. [ll]. We define a set of one- 
dimensional heights which reside on the vertices of the 
square lattice dual to the plaquettes without crossing in- 
teractions. Next, we assign an orientation (clockwise or 
anticlockwise) to each tetrahedron, so that neighbouring 
tetrahedra have opposite orientations. Since the lattice 
dual to the tetrahedra is bipartite, this can be done con- 
sistently The rule for the height differences is as follows. 
If going from one site of the height lattice to another one 
passes over an up (down) spin, one increases (decreases) 
the height by one provided the spin was crossed in the 
direction given by the orientation of the tetrahedra it 
belongs to. In the opposite direction, one decreases (in- 
creases) the height by one. This generates a consistent 
assignment of the heights since going around a unit cell 
of the height lattice generates zero height difference by 
virtue of the ground-state two-up two-down condition, 
and because each closed path on the height lattice can 
be decomposed into a combination of such elementary 
loops O 

We next consider the action of a longitudianl field of 
strength h in the absence of a transverse field. For h < 
2 J, all ground states remain degenerate since they have 
zero net moment. At h — 2J, the applied field is strong 
enough to surmount the exchange field and it generates 
a spin-flop transition to a manifold of states with three 
spins up and one down in each tetrahedron. These states 



continue to have an extensive entropy but one which is 
reduced compared to the low-field value. In addition, the 
classical correlations in this regime, to which we allocate 
the name IM ('intermediate'-field regime), continue to 
be critical. This result follows from another mapping of 
those states onto a dimer model which in turn generates 
a height model. 

This mapping is obtained as follows. Consider a dimer 
model on the square lattice dual to the tetrahedra. For 
each spin pointing down, place a dimer, centred on this 
spin, with its ends located at the points of the dual 
lattice denoting the centres of the tetrahedra the down 
spin belongs to. Since each tetrahedron has exactly one 
down spin, and since each spin is shared by two tetrahe- 
dra, each classical IM ground state generates a hardcore 
dimer covering of the (dual) square lattice and vice versa. 
Such a square-lattice dimer model can be mapped onto 
a height model, from which the criticality of the correla- 
tions in this regime follow. In the IM-phase, the entropy 
per spin is a quarter of the value found in the Villain 
model, as now there are four spins per dimer rather than 
one: S/k B = 0.073. 

We emphasize that the field leaves the entropy un- 
changed for a finite range of h and then reduces to a 
lower value without eliminating it completely, while in- 
ducing a transition between two critical states; this again 
persists over a finite range of fields. At h = 6J, there fi- 
nally occurs a transition to the fully polarised state. 

Next, consider tilting the field such that Y -C h; this 
enables us to generate a perturbation theory controlled 
by the small parameter Y/h. The transverse field induces 
matrix elements between the states corresponding to the 
classical ground state configurations. Both in the low- 
field and the IM regime, the degeneracy is not lifted until 
fourth order in perturbation theory. The reason is that 
connecting two ground states requires flipping a closed 
loop of alternating spins which passes through an even 
number of sites of any tetrahedron. The shortest such 
loop has length four and is thus generated at fourth order 
in perturbation theory. The lower-order terms induce 
only a diagonal shift in the energies but this shift is the 
same for all states. 

The short flippable loop, depicted in Fig. [ll], has differ- 
ent interpretations in the two phases. For the low-field 
phase, it corresponds to changing the height of a pla- 
quette whose four neighbouring heights are equal. In a 
manner analogous to the ground-state selection on the 
triangular lattice, this leads to the selection of the flat 
state in height language. In spin language, the flat state 
is a Neel state on the square lattice underlying the square- 
lattice with crossings. 

In the IM phase, the quadruple spin-flip again leads 
to the RK model at v = and t cx Y 4 /J 3 (see Eq. [O]). 
This move connects different ground states since the to- 
tal magnetisation of each tetrahedron remains unchanged 
and the effect of flipping the four spins is to generate the 
familiar dimer plaquette move I 1 <-> 

Carrying over the results from the fully frustrated 
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square lattice, we expect the system to order into a flat 
(height) phase which corresponds to a columnar dimer 
phase. In spin language, this phase differs from the Neel 
low-field phase, from which it can be obtained by flipping 
half of the down spins, e.g. those which are located on 
the top-left hand corner of one sublattice of tetrahedra. 

The complete phase diagram for this magnet in the 
h — r plane can be constructed from this. The simplest 
phase diagram incorpoating our results is displayed in 
Fig.Pl 



VIII. THE KAGOME LATTICE 

In the kagome Ising antiferromagnet, depicted in 
Fig. ||, the nearest-neighbor exchange couplings are uni- 
form and-,antiferromagnetic. The ground-state entropy 
is finite,E3 Skag — 0.502, and is more than half of the 
maximum paramagnetic value (S para = In 2) in con- 
trast to the triangular case (Stri = 0.323). Further- 
more Skag is close in .value to that obtained by the 
Pauling approximation,^^ S Pauling = In 2 + | In | w 
0.501, where the triangles are considered independently; 
this suggests that spin-spin correlations in the kagome 
IAFM are extremely weak. This is indeed the case, 
and the model remains classically disordered at all 
temperatures .E3 

An important feature or the kagome lattice is that the 
frustrated units (triangles) are arranged to share sites 
rather than bonds. This precludes the mapping of the 
kagome ground states to a hardcore dimer model on the 
dual lattice (known as the diced lattice); this also pre- 
cludes a mapping to a height model. It is thus an ex- 
cellent candidate for a disordered quantum magnet. The 
physics of this model has been discussed in moderate de- 
tail in Ref. |l6|; here, we fill in the missing detail, re- 
peating some material telegraphically for coherence and 
convenience. 

Application of a longitudinal field, \h\ < 4J, to this 
magnet leads to the development of a ferromagnetic mo- 
ment coexisting with a reduced but extensive entropy 
and critical spin correlations. The result follo ws in a way 
formally similar to the T = phases in Sect. VII. Each 
triangle has two up and one down spin, and denoting 
each down spin by a dimer with endpoints in the centres 
of the triangles that share it, we obtain a bijective map- 
ping between the ground states and the hardcore dimer 
covering of the hexagonal lattice; see Fig. [l^. This im- 
plies, by virtue of the associated height mapping, critical 
correlations. The entropy per spin, while still nonzero, is 
reduced to the value of S = Strife = 0.108; the differ- 
ence to the triangular lattice arises from the fact that the 
number of spins per dimer is different on account of the 
inequivalent mappings used to arrive at the dimer model. 

We emphasize that this result is rather unusual. Start- 
ing from a disordered magnet, we obtain a critical state 
with nonzero moment and extensive entropy upon appli- 



cation of an infinitesimal field; these properties, including 
the criticality, persist for a finite range of field strengths. 

Next, consider applying a small transverse field, T <C 
\h\, in addition to the longitudinal one (0 < \h\ < 4J). 
The resulting perturbation theory, controlled by the pa- 
rameter |r//i|, has the following structure. The ground- 
state condition of the exchange term of the Hamiltonian 
imposes the restriction of having either one or two down 
spins ('dimers') per triangle. Starting from a configu- 
ration which is also a ground state of the longitudinal 
field part of the Hamiltonian, we are thus allowed to add 
dimers violating the hardcore condition, as long as three 
never meet in one site. To connect different ground states 
requires relocating three dimers (denoted by crosses in 
Fig. |l^), since the shortest closed loop of bonds of the 
hexagonal lattice has length six. The lowest-order off- 
diagonal matrix elements thus arise to 0((T/h) e ), and 
are precisely those desc ribed by the quantum dimer res- 
onance term in Eq. 5.1. This being the unique shortest 
closed loop, all other terms up to and including sixth 
order are diagonal. These terms correspond to putting 
down and then removing up to three dimers. Due to the 
local structure of the dimer states, the number of such op- 
erations, and the concomitant energy denominators, are 
found to be the same for all ground states. Hence, the 
diagonal energy shift is uniform and does not generate a 
lifting of the classical degeneracy. 

The ordering patte rn i s therefore determined by the 
hexagonal QDM (Eq. |t]) at v = and t cx T(T/h) 5 . In 
dimer language, it is the one depicted in Fig. [7]. 

Next, consider the kagome IAFM in a transverse field 
r with h = 0. Following our previously described strat- 
egy, we look for a symmetry-breaking pattern within a 
Landau-Ginzburg- Wilson analysis; However, this mean- 
field treatment predicts an infinite number of zero modes, 
corresponding to the .simultaneous softening of an entire 
branch of excitations £3 For Ising spins, high-temperature 
series expansion studies of the kagome IAFM indicate 
that thermal, fluctuations fail to select a wavevector to 
any order .Ej EI 

The variational maximally fiippable configurations 
turn out to be the maximally polarised dimer configu- 
ration defined above for h > T = 0, of which there is an 
exponentially large number. The resulting hierarchical 
states are those in which each triangle has a spin point- 
ing up, one pointing down and one pointing along the 
transverse field. Evidently, these states map onto the 
three state Potts model on the kagome lattice, which is 
also known to have a nonzero entropy S. 

All these arguments portray the kagome TFIM as a 
system extremely reluctant to order. We have checked 
this explicitly by quantum Monte Carlo simulations and 
found that the TFIM mirrors the classical model, with 
correlations somewhat enhanced compared to the clas- 
sical ones, but still rapidly exponentially decaying (see 
Fig. 3 of Ref. |6|). Here, we supplement these data by dis- 
playing, in Fig. [l4|, that the simulated correlation func- 
tions have settled down with respect to both quantum 
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temperature and discretisation error. Note that in both 
cases the correlations are extremely small in magnitude 
below the first few neighbours, in marked contrast to the 
situation in models known to order; for example, in the 
triangular IAFM the saturated correlation function re- 
mains above 0.5 in these units at the largest distances. 

In Fig. [l5|, we also show that, for a small longitudinal 
field applied in addition to the transverse one, the rapid 
decay of the correlation functions remains unaffected, the 
main effect being the appearance of a net moment vis- 
ible in the correlations at large distances. The size of 
this moment is close to linear in h, as one expects for a 
quantum paramagnet. At short distances, we observe the 
emergence of weak correlations reminiscent of the dimer 
crystal described above. We therefore conjecture that 
with further field-tilting results in a continuous quantum 
phase transition to the ordered dimer phase. We have 
not been able to confirm this numerically as our simu- 
lations fail to equilibrate before the critical value of h is 
reached. 

In Fig. |l^ we display the simplest phase diagram for 
the kagome IAFM in longitudinal and transverse fields 
consistent with the discussion here, noting that details 
associated with the tilted-field phase line remains a topic 
for future study. We close this section by restating its 
main result, namely the fact that the kagome TFIM is a 
quantum disordered magnet. 

IX. THE SAWTOOTH CHAIN 

An extreme and amusing example of a quantum spin 
liquid is provided by the sawtooth chain (Fig. |l6|). It 
is the ultimate cooperative paramagnet: the locations of 
the frustrated bonds on each triangle are entirely inde- 
pendent, providing a ground-state entropy S = (ln3)/2 
per spin. This is a result of the absence of closed loops of 
triangles, and the fact that, as in the case of the kagome 
lattice, they are arranged to share corners. 

This chain has been studied in detail by Priour et 
al.Cj using a high-order series expansion in J/T, which 
compared favourably with exact diagonalisations. They 
found no phase transition at any value of J/T, implying 
that the chain is in a quantum paramagnetic state. 

It has been suggested that the ground-state topology 
plays an important pale in determining the ordering prop- 
erties of a magnet .Ea It therefore may be of interest to 
note that this chain has a completely connected ground- 
state manifold. The proof proceeds by explicit construc- 
tion of a path from any ground state configuration to a 
reference ground state configuration, pictured in Fig. |l^. 
First, one picks any down spin on a bottom row (if there 
are none, any one can be flipped). The remaining two 
spins of the triangle can be made to agree with the refer- 
ence state by one of the three operations depicted below 
the reference configuration; this is then repeated for the 
neighbouring triangle. In the case of the operation pic- 



tured on the left, nothing has to be done. A question 
mark means that the spin can have either orientation 
and can thus be made to agree with the reference orien- 
tation immediately. Once this has been done, the spins 
enclosed by the ellipse can be oriented appropriately. 

X. THE FULLY FRUSTRATED HEXAGONAL 
LATTICE 

A. The classical model 

The fully frustrated hexagonal Ising magnet is ob- 
tained from the ferromagnetic Ising model on that lattice 
by changing the sign of one interaction in each hexagon, 
as depicted in Fig. [l7|. Since the lattice is bond-sharing, 
the classical ground states can be again be represented 
by hardcore dimer coverings of the dual lattice, which 
is the triangular lattice. Unlike the triangular lattice, 
the fully frustrated hexagonal lattice does not admit a 
height model of the types encountered for other lattices 
we discuss. Although it is possible to give a rule to as- 
sign heights to the sites of the hexagonal plaquette, these 
rules cannot be consistently given for all plaquettes of the 
full lattice since they involve assigning opposite orienta- 
tions to neighbouring plaquettes, which is impossible as 
the dual triangular lattice is not bipartite. 

It is found that the classical correlations of the hexag- 
onal FFIM are disordered rather than critical, and the 
ground-state entraoy of the magnet has been evaluated 
to give S = 0.214E3 

B. Action of the transverse field and flippability 
analysis 

Let us now consider the dynamics induced by the trans- 
verse field. Note that the odd coordination of the lattice 
precludes any spin from being flippable individually since 
it cannot have an equal number of frustrated and satis- 
fied bonds. The shortest allowed dimer move consists of 
moving two occupied dimers and implies flipping a pair 
of neighbouring spins. Not absolutely all local moves 
within the ground-state manifold, however, can be gen- 
erated with these pair spin flips. In fact, one particular 
('staggered') configuration exists which does not allow 
any two spin flips but which instead permits the four 
dimer move as depicted in Fig. [D]. 

For infinitesimal T, the ground-state degeneracy is thus 
lifted to second order in T/J, and the relevant spin flips 
are the pair flips. To identify the maximally flippable 
configurations, we note that each dimer can be part of 
at most two flippable pairs. Since the total number of 
dimers is fixed, the maximally flippable configurations 
are those in which each dimer belongs to two pairs. 

All maximally flippable configuration can be obtained 
by carrying out any number of either operation A or B 
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(right panel of Fig. [Tt]) on a particular, maximally flip- 
pable ('columnar') configuration, in addition the sym- 
metry operations of global rotations and global inversion. 
These operations involve exchanging empty and occupied 
dimers on an infinite alternating sequence, along a string 
(A) or a sawtooth (B). They generate, for a system con- 
taining L 2 spins, a number of configurations exponential 
in L rather than L 2 . The configurations generated by op- 
eration A all incorporate long-range order in one special 
direction. 



C. Large T analysis 

Next, we carry out the large T analysis for this lat- 
tice. We first identify the soft modes, then construct 
the Landau-Ginzburg- Wilson action which we minmise 
in order to obtain spin configurations and correlations 
functions. 

Since the hexagonal FFIM has four sublattices (la- 
belled as in Fig. [IT]), the Fourier transform of the in- 
teraction matrix is of size 4x4, and is given by (omitting 
an overall factor of J/2: 



V 



1 - e lq * 
- e~ l i* 

, e -iqy 



. e l Qv 











(10.1) 



The square of the eigenvalues, A, are 



3 ± y 6 + 2 cos^^) - 2 cos(q x + q y ) + 2 cos(q x - q y ) 

The four extremal eigenvalues of interest are V6, which 
occur at wavevectors ±(w/6, n/2), ±(57r/6, ir/2). The 
corresponding eigenvectors occur in complex conjugate 
pairs: 
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where F = 2 sin(57r/12). 

In order to determine the terms in the LGW Hamilto- 
nian, one determines how these modes transform among 
themselves under the symmetry operations (translations 
(T x ,T y ), reflection (Ref) and rotation (Rot)) of the un- 
derlying lattice. Each of these symmetry operations 
comes with a gauge transformation since the unit cell 
of the lattice and the unit cell of the interactions are un- 
equal. These are shown in Fig. |1| The transformation 
matrices for the amplitudes of these modes under the 
abovementioned symmetry operations are thus: 
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To find the eventual action, one has to determine the 
terms at each order which remain invariant under these 
transformations. One obtains: 



C={r + q 2 )^ 2 a + ^ 2 ) + u^ 2 a 
+u 6 ^ 2 a + ^ 2 b ) 3 + 
+v 6 (ip a if)l cos(6> a - 56» 6 ) 
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Here we have transformed the four mode amplitudes 
{a.i\i — 1...4} into complex numbers il) a ,b exp(i# a &) with 
ip a exp(z# a ) = ai + ias and ipb exp(z6'b) =02+ ia^. 

The terms which survive are 0(4)-symmetric up to 
sixth order, where in addition a symmetry-breaking term 
appears. To minimise this action, we use the 0(4) sym- 
metric terms to determine ip 2 + ip 2 , the radius R of the 
0(4)-sphere, which sets the overall amplitude of the spin 
pattern. The symmetry-breaking wg-term then selects 
48 points on this sphere. One finds that ip a — ip + and 
ip a = ip+/F; 8 a , b = (±7r/6) for v 6 < (>) 0. The other 
47 solutions are generated by the operations ip a <-> V^fc 
and {8 a -> 6 a + tt/12, 9 b -> 9 b + 5tt/12}. The 48 solu- 
tions for a given sign of vq are related by simple symme- 
try operations. Note the large unit cell of the ordering 
pattern, which contains 48 spins. This can in part be 
attributed to the non-uniform pattern of bonds. In the 
gauge-invariant dimer language, the unit cell could be 
smaller; an example of this is given further down. 

The correlation functions obtained from the above ex- 
pressions are for soft spins so that they will certainly 
not be quantitatively found in the real system at low 
temperatures. However, their qualitative features, if 
the ordering pattern is correctly predicted, should sur- 
vive. These are peaks at the appropriate wavevectors in 
Fourier space; due to the factors of 1/F for the ampli- 
tudes ij) a ,bi these peaks are not expected to have equal 
heights for correlation functions for spins on the same 
sublattice. It turns out that the calculated correlation 
functions when, averaged over all 48 minima, are the 
same as those averaged over the entire 0(4) sphere. 

We have looked for this ordering pattern by quan- 
tum Monte Carlo. There are several features complicat- 
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ing this search. One arises from the quantum dynamics 
which takes the form of double spin flips. Luckily, one is 
saved here by the fact that the cluster algorithm can be 
generalised to this case. The most simple-minded exten- 
sion to double spin flips would have been to generate two 
neighbouring clusters (rods in the imaginary time direc- 
tion) separately - by design, the cluster algorithm would 
cancel the problematic Boltzmann factor exp(— K T /2) in 
each of the two clusters separately. However, clusters 
of unequal height, which become overwhelmingly proba- 
bly for large K T , cannot be flipped because permissible 
moves must be double spin flips in any and all planes. 
The way around this is to construct a cluster as a pair 
of rods. This works because the cluster can either ter- 
minate when encountering zero, one or two domain walls 
in the time direction. The case of zero and two can be 
taken care of by working with an effectively doubled cou- 
pling in the time direction; the case of one domain wall 
takes care of itself as there is no Boltzmann factor to be 
canceled in the first place. This prescription can there- 
fore not be generalised to other multiple-spin flips in a 
straightforward manner. 

Another problem is the very large unit cell of the or- 
dering pattern, which is 48 spins but only fits into pe- 
riodic boundary conditions for system sizes multiples of 
192 spins. We have thus not attempted to do a complete 
finite-size scaling study and have contented ourselves 
with displaying the presence of peaks at the wavevectors 
predicted above. 

In Fig. [l^ we have displayed the correlation function 
for spins on sublattice 4 along the x and y directions (also 
see Fig. 17). A peak at the expected location is clearly 
visible. This suggests that there is at least significant 
short-range order of the predicted nature present in this 
regime. We caution, however, that we have not been able 
to do a complete numerical analysis of this problem that 
would unambiguously establish this ordering pattern. 

To rationalise the above results, we finally present a 
dimer pattern which corresponds to a spin pattern of the 
correct size unit cell (Fig. ^(J) .cJ We have not found a way 
of relating the spin-spin correlations precisely to those of 
the solutions of the LGW theory found above, but one 
can see rather nicely how the periodicity of twelve and 
four in the x and y directions comes about from a smaller 
size unit cell in the gauge invariant description. The 
state pictured is favourable because each spin fluctuates 
in a pair gaining energy from the transverse field, and 
because there are many equivalent such configurations, 
as explained in the caption. 



of theoretical interest as the simplest settings in which 
quantum dynamics interacts with classical frustration. 
They can represent effective theories of systems, such as 
short ranged RVB magnets, where the low energy dy- 
namics contains a frustrated Ising degree of freedom and 
are likely of experimental interest when suitable materi- 
als are probed by the application of a transverse magnetic 
field. 

We have been able to make considerable progress in 
elucidating the phase structure of these systems which 
sheds much light on the interplay between the structure 
of the macroscopic degeneracy and the quantum dynam- 
ics. Our results include instances of "order by disor- 
der" and of "disorder by disorder" (i.e. Ising spin liq- 
uids). In arriving at these we have used two system- 
atic approaches: a variational approach that builds on 
the local entropy of different classical configurations that 
is the Ising analog of semiclassical analyses of magnets 
with continuous symmetry, and an LGW analysis that 
attempts to guess at the large order structure of an ex- 
pansion about the large transverse field paramagnet. We 
have also made use of quantum Monte Carlo simulations 
as well as of mappings to height and dimer models with 
the latter allowing a connection to RVB physics. 

There appear to be several directions that can be pur- 
sued further. Along the lines of the questions addressed 
in this paper, the analysis of the kagome system in tilted 
fields and a definitive analysis of the hexagonal lattice 
problem are needed. Further, he structure and energet- 
ics of the low energy excitations needs to be worked out 
to have an understanding of the dynamics at low tem- 
peratures. Beyond this it would be interesting to include 
the dynamics of the transverse exchange and see what 
happens both in the approach to the Heisenberg limit 
and when a transverse field is also present. We expect 
such work to be both fruitful and instructive. 
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XI. SUMMARY 



In this paper we have described our analysis of a num- 
ber of frustrated Ising systems with a quantum dynam- 
ics introduced by a magnetic field applied transverse to 
the Ising axis. We have argued that these models are 
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FIG. 1. Two frustrated plaquettes with Ising spins. 
In the canonical example of the antiferromagnetic triangle, 
anti-aligning two spins leaves the direction of the third unde- 
termined. Similarly, a square can be frustrated by choosing 
an odd number of bonds to be antiferromagnetic. For such 
mixed-bond models, we represent (anti)ferromagnetic bonds 
by (dashed) solid lines. 







FIG. 2. The lattices on which fully frustrated trans- 
verse field Ising models are discussed in this paper. Clock- 
wise from noon: kagome, hexagonal, pentagonal, saw- 
tooth chain, three-leg ladder, square lattice with crossings 
("two-dimensional pyrochlore"), square and triangular. For 
normally unfrustrated lattices, frustration is introduced by 
choosing an odd number of bonds in each plaquette to be 
antiferromagnetic (see Fig. |l|). 



FIG. 3. The fully frustrated three-leg ladder. Solid lines 
correspond to ferromagnetic bonds, dashed lines to antiferro- 
magnetic ones. Dimer representations of the columnar (right) 
and one of the two staggered (left) configurations. 
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FIG. 4. Ground-state energy (proportional to (S x )), en- 
ergy gap and {S z ) for the three- leg lader in the limit r — » + 
from exact diagonalisation. 
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FIG. 5. The phase diagram of the Rokhsar-Kivelson quan- 
tum dimer model for the fully frustrated Ising magnet on the 
three-leg ladder. The diagonal correlations at v = t are those 
of the classical dimer model. The infinitesimal transverse field 
problem, T — + , maps to the point v — and thus corre- 
sponds to a finite jump of v. 
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FIG. 6. The pentagonal lattice with its maximally flippable 
configuration. The three spins in the unit cell are labelled by 
(0,1,2), and the x and {/-lattice translation vectors are given 
by the arrows. The leading-order perturbation theory reduces 
to studying the central horizontal ladder, which is equivalent 
to the three-leg ladder of Fig. H 




FIG. 7. The maximally flippable spin state (arrows) and 
its corresponding height and dimer configurations on the tri- 
angular lattice. The result of a single spin flip for the dimer 
and height configurations is also shown. See the text for de- 
tails. 
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FIG. 8. The phase diagram of the triangular TFIM. 
pattern of the ordered phase is depicted in Fig. (?) 
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FIG. 9. The absolute value of the Fourier transform 
of the spin-spin correlation function along the ^-direction, 
(S(x)S(O)) , (arbitrary units, at low temperature), as a func- 
tion of q x ranging from to n. Note the growing and narrow- 
ing Bragg peak. 
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FIG. 10. The correlation matrix, (MiMj) = (MjMi) of 
sublattice magnetisations, Mi, for a three sublattice state, 
as a function of the quantum temperature for a system with 
729 spins. The sublattices are labelled by the size of their 
magnetisation so that M x > M 2 > M 3 . (M|) ~ implies 
that one sublattice has zero root-mean square magnetisation. 
(Ml ) ~ (M|) ~ - (M1M3) implies the other two sublattices 
have equal and opposite magnetisation. 
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FIG. 11. The square lattice with crossings, or 

two-dimensional pyrochlore. A spin and the resulting height 
configuration are shown. Other features are explained in the 
text. 
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FIG. 12. The phase diagram in the h — F plane for the 
square lattice with crossings. The precise extent of either 
phase away from F/h <C 1 is unknown. 
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FIG. 13. Phase diagram for the kagome Ising antiferro- 
magnet in a field. Inset: Mapping of the kagome IAFM in 
a longitudinal field onto the hexagonal lattice dimer model. 
The down spins are marked by dimers. The non-trivial move 
to sixth order in perturbation theory corresponds to flipping 
all the spins marked by crosses. 
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FIG. 14. The angularly averaged spin-spin correlation 
function for different quantum temperatures and discretisa- 
tions, A = exp(K T /2). Note the rapid decay and the enlarged 
scale. The quantum curves lie almost on top of one another. 
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FIG. 15. The angularly averaged spin-spin correlation 
function for the kagome TFIM in a longitudinal field. The 
correlations corresponding the dimer crystal (Fig. [l3]), scaled 
down by a factor of 10, are also shown. Note that the 
high-field curve starts moving towards the crystal correlations 
at small distances. Inset: The magnetisation per spin (ordi- 
nate) is linear in field {h/F, abscissa). 





FIG. 16. The sawtooth chain with the reference 

ground state and the three operations needed for proving 
ground-state connectedness (see text). 
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FIG. 17. The fully frustrated hexagonal lattice. The ba- 
sis spins of the lattice are numbered 0-3 and the rectangu- 
lar lattice translation vector in the x (y) direction is given 
by the horizontal (vertical) arrow. The pair spin flip and a 
non-trivial multiple spin flip are indicated with their corre- 
sponding dimer moves. Maximally flippable states are gener- 
ated from the columnar seed state by interchanging occupied 
and empty dimers along horizontal lines (A) or in a sawtooth 
pattern (B). 




FIG. 18. The symmetry transformations and concomitant 
gauge transformations for the Landau- Ginzburg- Wilson ac- 
tion. Translation in the ^-direction (arrow in left panel) 
require no gauge tranformation. Reflections about the 
dot-dashed line go along with inversion of the spins marked 
by crosses; translations by the other hexagonal lattice vector 
(arrow in the right panel) require flipping the circled spins, 
and rotations by n/3 (circle with arrow) involve flipping the 
spins denoted by ticks. 




FIG. 19. The Fourier transform of the absolute value of 
the autocorrelation function of a spin on sublattice 4 in both 
x and y directions for a system of 4608 spins. Results for the 
classical and the quantum (F = + ) cases are shown, with 
a quantum temperature Tq/T = 0.028 nominally, but the 
discretisation error is rather large: A = 4.5 only. 




FIG. 20. Cartoon of a dimer pattern with a maximal num- 
ber of independently flippable pairs of spins. To obtain the 
corresponding spin pattern, one has to choose a gauge, i.e. one 
has to make a choice of which bonds are antiferromagnetic. 
Fig. [I?] shows one such choice. For this choice, the spin pat- 
tern has a 48 site unit cell (two plaquettes on the left), whereas 
the gauge invariant dimer description only has 24 sites in a 
unit cell (any single plaquette). A rhombus stands for a pair 
of dimers (frustrated bonds) resonating between two configu- 
rations, so that the two spins enclosed by it can be thought of 
as flipping together to take advantage of the transverse field. 
This pattern cannot be derived from a uniform maximally 
flippable state but has extensive 'configurational' entropy, as 
there are two distinct ways (shown in the left and right pairs 
of plaquettes, respectively) of pairing up the spins in each 
plaquette. 
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